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Abstract 



We describe a Bethe ansatz based method to derive, starting from a mul- 
tiple integral representation, the long-distance asymptotic behavior at finite 
temperature of the density-density correlation function in the interacting one- 
dimensional Bose gas. We compute the correlation lengths in terms of so- 
lutions of non-linear integral equations of the thermodynamic Bethe ansatz 
type. Finally, we establish a connection between the results obtained in our 
approach with the correlation lengths stemming from the quantum transfer 
matrix method. 
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1 Introduction 



This article is devoted to the study the long-distance asymptotic behavior of the correlation 
functions in the one-dimensional Bose-gas at finite temperature. The method of analysis we 
apply builds on the method developed in the work [IJ. There, the long-distance asymptotic 
behavior of the zero temperature correlation functions of massless quantum integrable models 
was derived in the framework of the algebraic Bethe ansatz. This setting allows one to present 
the correlation functions at zero temperature as series of multiple integrals |2] |3l IU [5l |6]. In 
the long-distance regime, each of the multiple integrals appearing in one of these series can be 
evaluated asymptotically. The resulting asymptotic series can then be summed up [Ij. 

In [1], we have focused on the example of the generating function for correlation functions of 
the third components of local spin in the XX Z spin-1/2 Heisenberg chain. We have mentioned, 
however, that the method can be applied to other integrable models as well. In particular we 
showed how one can use the same approach for the model of one-dimensional bosons described 
by the quantum nonlinear Schrodinger equation (QNLS model) using [6j. 

We would like to stress that the algebraic Bethe ansatz approach of [1] was fit for the 
analysis of the asymptotic behavior of correlation functions at zero temperature. However, in 
some models such as the QNLS one, it can also be applied for the evaluation of the long-distance 
asymptotic behavior of correlation functions at finite temperature. It is this problem that we 
consider in the present article. 

The Hamiltonian of the QNLS model is given by 



coupling constant and h the chemical potential. We focus on the case of the repulsive regime (i.e. 
c > 0) and consider the model on a finite interval [0, L] subject to periodic boundary conditions. 
We will take the thermodynamic limit L — t- oo later on. 

This model, for generic c, was first introduced and solved by Lieb and Liniger in 18]. It can 
be considered as a generalization of the model of impenetrable bosons considered by Girardeau 
[9]. The spectrum of this model can be obtained by the algebraic Bethe ansatz |10| \TT\ 112) . 
The thermodynamics of the QNLS model were first studied in [13J and this analysis was made 
rigorous in [MJ. It was shown there that the state of thermal equilibrium is described by a 
non-linear integral equation for the excitation energy e(A) 



L 




(1.1) 







Here ^ and are Bose- fields subject to canonical, equal-time commutation relations, c is a 



e{X) = X^ -h 
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where T is the temperature and 
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Below we refer to equation ( 1.2 ) as the Yang- Yang equation. The state of the thermal equilibrium 
is given by a Dirac sea filled with a certain density of particles Pp{\) and an associated density 
of holes The total density is denoted pt(A) = Pp{X) + PhW- The Fermi weight i?(A), as 

usual, is defined as the ratio of the density of particles to the total density, namely, 

^-^')-^(A) = (l + e^)'\ (1.4) 



and the functions Pp{\) and pt(A) are related to each other by the integral equation 

If 1 

The temperature dependent correlation functions are defined in a standard way 

tr (0e-^/^) j:{n\0\n)e-^/^ 



Here |r2) are eigenstates of the Hamiltonian ( |1.1[ ) and E are their eigenvalues. The sums in (1.6) 
are taken over the complete set of the eigenstates 

Exact representations for the temperature dependent correlation functions of the QNLS 
model where obtained for the impenetrable Bose gas (c = oo) in |15| [THl [T7] by using the 
free fermion structure of the model. The long-distance asymptotic behavior of the two-point 
functions in the impenetrable Bose gas has been derived in |18| [T9l |20] with the use of Riemann- 
Hilbert problem method. The case of general coupling constant < c < -|-oo was considered in 
|21| by the algebraic Bethe ansatz. There a series of multiple integrals was obtained for thermal 
density-density correlation function of the QNLS model. The method of dual fields was applied 
in |16| W2\ [23l for the derivation of determinant-type representations for various correlation 
functions. Those representations were used in |251 [26] for the asymptotic analysis. One should 
also mention the method of the asymptotic analysis based on the functional integral approach 
|27| . Yet another approach to estimate these long-distance asymptotic behavior of correlation 
functions at finite temperature is provided by the conformal field theory )28l |29| l30l [3T] and 
references therein. These last two methods are, however, restricted to the low-temperature 
regime only. 

Recently, in the work [32], the quantum transfer matrix (QTM) approach was applied to the 
QNLS model. We recall that, originally, this method was developed for quantum spin chains 
(see the nice review [33] and references therein), where one can construct the QTM Tg explicitly. 
The diagonalization of the QTM gives access to its leading eigenvalue (the logarithm of which 
corresponds to the free energy) as well as to the subdominant ones (which in their turn give 
access to the correlation lengths at T > 0). In the infinite size limit, these eigenvalues are 
expressed as weighted integrals involving a counting function which satisfies a thermodynamic 
Bethe ansatz (TBA) non-linear integral equation (see also |34]). To the best of our knowledge, 
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the analog of Tq in the QNLS model is not known nowadays. Therefore, in [32| . the last model 
was treated as some special, continuous limit of the XXZ spin chain. It was shown that, in 
this limit, the non- linear integral equation describing the maximal eigenvalue of the QTM goes 
to the Yang- Yang equation. This observation allowed the authors of [32] to obtain a multiple 
integral representations for the temperature correlation functions in the QNLS model. 

The first steps of our method are close to the ones adopted in |21]. Namely, we build on the 



arguments given in |35j (see also \21\ I12|): in the thermodynamic limit, the representation (1.6) 
can be replaced, at least for local or quasi-local operators O, by a single expectation value 

where \^t) is any one of the eigenstates of the Hamiltonian corresponding to the thermal equi- 
librium. This representation constitutes the starting point for our calculations. The matter is 
that the algebraic Bethe ansatz provides multiple integral representations for the expectation 
values of a wide class of operators with respect to an arbitrary eigenstate of the Hamiltonian. 
In particular, one can obtain such representations for the expectation values with respect to the 
state of thermal equilibrium. In this way, we recast {0)t as a series of multiple integrals. Once 
this is done, it remains to evaluate the asymptotic behavior (in the distance) of the obtained 
integrals. It is remarkable that our results are given in terms of solutions to TBA non-linear 
integral equations. These can be understood as the limiting (in the sense of |32) ) equations 
describing the sub-leading eigenvalues of the QTM. Thus, our method provides a link between 



the method based on the representation (1.7) and the QTM approach. In particular, we obtain 
an extension of the QTM-based results |37| |38] for the correlation lengths in the XXZ spin chain 
here for the QNLS model. By applying the scaling proposed in [32j one can actually map the 
correlation lengths obtained for the XXZ spin chain to the ones obtained by our asymptotic 
analysis. 

We have already mentioned that the representation for temperature correlation functions of 
QNLS model as a series of multiple integrals were obtained also in |21] . Later it was analyzed in 
[39\ HOI HI] . However the authors of those works restricted themselves to the analysis of the first 
few terms of the series. This led them to the wrong conjecture for the long-distance asymptotic 
behavior of the correlation functions. We would like to stress that the series of multiple integrals 
obtained in [2T] . as well as ours are not well ordered in respect to the x — )• oo limit. That is to 
say all the summands in both series contribute to the leading asymptotic behavior as well as to 
the corrections. Our method allows us to compute all such contributions and re-sum them. 

This article is organized as follows. In section [2] we state the problem to solve and the main 
results of this paper. In section [3] we derive a series of multiple integral representation for the 
generating function of the density-density correlation function at finite temperature. In section [4] 
we perform the asymptotic analysis of the individual multiple integrals appearing in the series. 
We re-sum these asymptotic expressions in section |5] by the use of Lagrange series. Then, in 
section [6] we obtain the long-distance asymptotic behavior of the density-density correlation 
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function. There, we discuss the leading term and connect the various correlation lengths we 
obtain with the quantum transfer matrix method. In appendix[A]we present a numerical analysis 
of the expressions for the correlation lengths that we have obtained. Technical aspects relative to 
the formulae for the amplitudes are gathered in appendix|B] We discuss certain integral identities 
in appendix |C] Finally, in appendix |D] we briefly review the results we need on continuous 
Lagrange series. 



2 The problem to solve and the main results 

To calculate the long-distance asymptotic behavior of the density-density correlation functions 
we consider its generating function 

were a is a complex numbei[^ This generating function is defined in terms of the number of 
particles on [O5 x] operator Qx 

X 

Qx = jjiz)dz, with j{x) = ^\x)^{x), (2.2) 



Then, the density-density correlation function {j{x)j{0))T is obtained from (2.1) by 

-1 92 Q2 



(2.3) 

Q = 



To describe the large-x asymptotic expansion of (e^'^*"®^)^, we first consider a set of functions 
Ui{X) satisfying the non-linear integral equations 

Ui{X) = X"" - - ^ j K{X-fi)\og(l + e-^) dX - iTj2[^{s+ - X) - 9{sj - X)], (2.4) 

TCP i=l 



where n = 0, 1, . . . , the kernel K{X — fi) is given by (1.3), ha = h-\- 27riaT^ and 

^(A) = i log ( , 0'W = K{X). (2.5) 

The parameters {sf} belong to the upper (resp. lower) half-plane and are the roots of the 
equation 

l + e-«.(^?)/r = 0. (2.6) 



^We draw the reader's attention to the fact that the combination 2Tvia was denoted by /3 in p]. We did not 
use this notation here so as to avoid a confusion with the inverse temperature, that is traditionally denoted by /3. 
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The subscript i in Ui serves as a way to label all the possible choices of n roots {sf}i. In 
the framework of our analysis, the roots sf arise as deformations of the poles rf of the Fermi 
weight (1.4). Namely, there exist functions sf{'y) of some parameter 7 such that sf{0) = rf 
and sf{l) = sf. Thus, one can associate every set {sf} with its pre-image {rf}. Every choice 
of the {rf}i uniquely yields the image {sf}i and the corresponding function Ui{X). Thus, the 
subscript i in the equation (2.4), in fact, enumerates different subsets of the Fermi weight poles 
{rf}i. Note that, for a given function Uj(A), there may exist other points Wp 7^ sf such that 
1 + e~"'('"p)/-^ = 0. Following the terminology of the QTM-based approach, we call such roots 
hole-type solutions. 

Equation (2.4) differs from the Yang- Yang equation (1.2) by a shift of the chemical potential 
and the term depending on the functions 9{sf — A). Due to (2.6), one can get rid of this term 
by changing the integration contour. Namely, let the contour d be a deformation of the real 
axis such that moving from M to d one only crosses the roots {sf}i while all other solutions Wp 
as well as all the poles of the Fermi weight rf are not crossed (see Fig. [T|. Then, equation (2.4) 
turns into 



Ui{X) = \^ - ho 



T 
2^ 



i^(A-/i) log 1 + e 



"i(M) 
T 



dA, 



(2.7) 



and has the form of the Yang- Yang equation up to the shift of the chemical potential and the 
change of the integration contour. Hence, different functions Ui{\) are enumerated by different 
contours C^: Ui{X) = u{X,Ci). It is interesting to notice that this method of recasting equation 
(2.4) into (2.7) is very similar to the one used in |42) to generate higher excited states within 
the TEA equations by analytic continuation. 




Figure 1: The poles of the Fermi weight rf are denoted o, the zeros sf by •, the hole- type 
solutions Wp hy x . The integration contour Cj bypasses the points s^ , S3 in the upper half-plane 
from above and bypasses the points s^, §2 in the lower half-plane from below avoiding the points 
rf and Wi. 
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The large x asymptotic expansion of ^e^'^*"®'")^ has the following form 

^g2™Q.^^^^g-xp,^[^^]^ x^oo, (2.8) 

i 

where 

Pr = -^ log ^ + " " dA , (2.9) 



i V 1 + e — / 

and the constant in x amplitudes B\ui\ are functionals of Ui{\). Their explicit but rather 



cumbersome form is given in section [6| It follows from (2.3) that the correlation lengths pi at 
a = drive the long-distance exponential behavior of the two-point function. Our numerical 
computations (c/ appendix [A|) together with the low-temperature expansions |47) confirm that 



given a = and any i, i.e. for any solution Ui to (2.7), K(pj) > and even = if Cj = M and 
5R(pi) > for all other contours Ci ^ M. 



The sum in (2.8) is taken with respect to all possible solutions Ui{\), or what is equivalent, 
with respect to all different contours Ci including Ci = M. In the last case, the function Ui{\) is 
equal to the excitation energy e(A) with the chemical potential shifted by 27riaT. It is precisely 
this term that gives the leading contribution to the asymptotic behavior of the density-density 
correlation function (see section [6]) . 



In (2.8), we did not indicate the upper limit of the summation as well as the possible cor- 
rections. We postpone the discussion of these topics to section [6] Here, we would like to draw 
the reader's attention to the obvious analogy of our results for the correlation lengths with the 
ones obtained by the QTM approach for the temperature correlation functions of spin chains 
|33| . We stress however that our method is completely different and that it is applied directly 
to the QNLS model without any limiting procedure from the XX Z Heisenberg chain like the 
one in 



3 Multiple integral representation 

The master equation [H |6] is a single multidimensional contour integral representation valid for 
a wide class of integrable models. It allows one to obtain various types of series of multiple 
integral representations for certain correlation functions. In the case of the one-dimensional 
Bose gas at finite temperature, a series of multiple integral representation was obtained in |32| 
through a scaling limit from the one of the XXZ chain. Here, we work with another one that we 
obtain from |6j by literarily repeating the steps given in |lj for the case of the XXZ Heisenberg 
spin chain. 

Prior to taking the thermodynamic limit, one can describe the eigenstates \Q) of the Hamil- 



tonian (1.1) by sets of real parameters A = {Ai, . . . , Atv}, A'' = 0, 1, ... , i.e. |0) = |r2(A)). These 
parameters \j are solutions to the Bethe equations |1H IT2] , and different sets A yield different 
eigenstates |r2(A)). 
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We remind that the master equation refers to a multidimensional contour integral representa- 
tion for the normalized expectation value of the e^'^*"®'" operator with respect to any eigenstate 
1 0(A)) of the Hamiltonian. Moreover, the expansion of this multidimensional contour integral 
representation into the aforementioned series of multiple integrals does not depend on the spe- 
cific choice of the state |il(A)). Hence, even when computing the thermal averages, one can still 
apply the results of [1]. Therefore, 

(Q(A)|e2™2.|5^(A)) _ 1 r^dzj ^ x- M 



MA)|J](A)) det^e^ n! J H 2vri ^ V M 



xdetf— ^ TT fV^^ > (3.1) 



where 



N 



2nLpiXj) = L + Y, K{\, - AJ, Q,k = Sjk - ^["^^.^^y . (3.2) 



a=l 



The functions Tn appearing in (3.1) are symmetric functions in the n variables {z} and in the 
n variables {A}. Their detailed description will be given later on. What only matters for the 
moment, is that these functions are holomorphic with respect to every variable in some open 
strip of fixed width a around the real axis. As for the discrete summations, each parameter /x^ 
runs independently through the entire set A. The integrals over Zj are taken with respect to 
a counterclockwise oriented, bounded and closed contour r(A) surrounding the set A in such a 
way that the only singularities of the integrand within the contour are the poles at zj G A. In 
particular, the contour r(A) lies inside of the strip of width a. 

In order to obtain the thermodynamic limit N,L ^ oo, N/L = D = const of the expectation 
value of the e^'^*"®"' operator, one starts by choosing a particular set A and then introduces the 
density of particles defined by Pp^{Xj) = limN,L^oo L{Xj+i — Xj). In the QNLS model, it can 
be shown that for sets A having a physical interpretation (the ground state, excited states of 
finite energy above the ground state, states of thermal equilibrium) this limit always exists. The 
density of particles allows one to replace the discrete sums over pk G A by integrals, 

hm y J;/(m)= f f{X)p,{X)dX. (3.3) 

Above, we have made the assumption that f{X)pp{X) is integrable on M. 

We now consider the case where the state |rj(A)) corresponds to any finite A^, L representative 



of the state of thermal equilibrium \0,)t- Then, due to (1.5), it follows that the thermodynamic 
limit of the function p{X) coincides with the total density pt{X). Respectively, the determinant 
of the matrix 0^^ goes to the Fredholm determinant 

lim det G,fc = det [/ - J-K'^'^ , (3.4) 
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where we used (1.4) and have introduced the kernel K^^\X, fi) = 'd{fi)K{X — fi). Finally, using 
(3.3) we can replace the discrete sums over the set A by integrals in every n*^ term of (3.1). 
Then agreeing upon 



we obtain 



E: 

ra=0 



(27ri)" 



»T = (e 



2TTiaQx ^ 



T ■ 



det 



(3.5) 



r({A}) 



{2-Kiy 



n 



ix(\j — Zj) 



^3 ~ ^3 



{z] ) n \Zk - Xj 



(3.6) 

, An avoiding any other 



where the contour r({A}) surrounds counterclockwise the variables Ai 
singularity of the integrand. Observe that apart from the replacement of the discrete sums by 
integrals we also have replaced the finite sum over n in (3.1) by the infinite series in (3.6). The 
question of convergence of this series will be discussed later, after the description of the functions 
J-n entering the representation (3.6) for ((e^'^*"®''))^ will be given. 

These functions depend on 2n variables Ai, . . . , An and zi, . . . ,Zn and read. 



n 



{A} 
{4 



{A} 
{4 



n^"(A 



I {A} 



(3.7) 



The explicit formulae for the functions Wn and Vn involve the set of auxiliary functions: 



VaA^^) = n 



a=l 



fi — Xa + ica 
fj- Za + ica' 



(T = 0,±, and -^^^(A) 



1 



X + ic X — ic 



Namely, 



{A} 

{z} 



The expression for Wn is more involved. It can be represented as 

y-,n{zk) 



Wn 



where 



(e 



{A} 



WK;n] ■ n 



k=l 



l)2det[I+^C7W]det[/+2^t/(2)] 



V_ 



VlAOl)) {V-AG2) - e2--VV;„(e2)) 



Above appear two Fredholm determinants of integral operators whose kernels are 



Ka{w - w') - Ko^iOi - w') 



(3.8) 



(3.9) 



(3.10) 



(3.11) 



(3.12) 



and 



Tj(2)( l\ T/ ( '\ Ka{w - w') - Ka{w - 62) 



(3.13) 
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These operators act on a counterclockwise oriented contour r(M) surrounding the real axis. By 
definition, the only singularities of U^^\w, w') inside of r(M) are the zeros of Vb;n(tL') (resp. the 



poles of Vb;n(u^'))- the following, when we will slightly deform the integration contour in (3.6), 
we will always keep this prescription for r(M), making it larger, if necessary. The parameters 6i 



and O2 in (3.11 )-(3.13) are arbitrary complex numbers lying inside of the contour r(M). It was 



proved in |T] that the overall combination appearing in the r.h.s. of (3.11) does not depend on 
a specific choice of these parameters. 

Due to the rather complicated form of the functions we are unable to provide a proof 



of the convergence of the series (3.6) in the case of general c. In fact, the situation is here 



quite analogous to the one occurring in the zero temperature case [1]. One can easily prove the 
convergence in the special case corresponding to the free fermion point c = 00. Indeed, then 



-,2ma 



1)" and the integrals over Zj can be taken explicitly. We obtain 



2-KiaQx 



E 

n=0 



2ma 



det 



sm 



j(Aj — Afc) 
7r(Aj - Afc) 



n 



(3.14) 



The series (3.14) is an expansion of the Fredholm determinant of the integral operator / + Vq 
with 



V^{X,^,) = {e 



2TTia 



1) 



sm 



7r(A- ^) 



(3.15) 



The general theory of Fredholm determinants ensures that the series (3.14) is absolutely con- 
vergent and defines an entire function of e^'^*". It particular, the expansion coefficients decay 
faster than exponentially. 



In the case c < 00, the series (3.6) cannot be reduced to a simple form as in (3.14). Nev- 



ertheless, taking into account that the QNLS model with a general coupling constant can be 
considered as a smooth deformation of the free fermions case, we shall assume in the following 



that the series (3.6) is absolutely convergent. 



In order to study the series (3.6) in the next sections, let us make two important remarks 



concerning the functions Wn, Vn and J^n present in (3.6). Although their precise expressions 



are rather cumbersome, it is easy to show that Wn and V„ (and hence Tn) ai'e holomorphic 
functions within some multi-dimensional strip (3: |9(2;j)| < a, |9(Aj)| < a, j = l,...,n. 
Though we can not determine explicitly the width a of this strip, it is clear that the latter 



is temperature independent, because the Fermi weight 'Q{X) (1.4) is the only function present 



in (3.6) that depends on T. In the case of positive chemical potential, the poles of the Fermi 
weight accumulate, in the T — )• limit, on certain points of M. Hence, there exists a crossover 
temperature Tq such that, for T < Tq, there are always poles of the Fermi weight that are located 
inside of the strip |9(A)| < a. The precise value of Tq depends on the value of a which, in its 
turn, is fixed by the analytic properties of J^n- For instance, in the free fermion point, one has 
Tq = -|-oo as J^n = (e^'^*" — l)"". We will not study this question further. Simply, when c < 00 
we shall limit ourselves to the regime T < Tq. In such a case, the nearest (to the real axis) 



singularity of the integrand in (3.6) always correspond to the poles of the Fermi weight. 
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The other important features of the functions Wn and V„ concerns their reduction properties, 
namely, 



V„ oj 



{A} 



\j=Zk 



{A} 



n-l 



{A} \ A, 

{z} \ Zk 



(3.16) 



In the following, we will see that these reduction properties (3.16) are crucial for the re- 



summation the series (3.6) in the asymptotic regime x — )• oo. 



4 Asymptotic behavior of multiple integrals 



The asymptotic analysis of the series ( |3.6[ ) can be done along the lines of the zero-temperature 
case [I]. The first step of that method consists in extracting the large-x asymptotic behavior 



of each 2n-fold integral appearing in the series (3.6). The main difference between the T = 



and T > situations is that, in the first case, the corrections to the leading terms have a 
power-law behavior in x, whereas, at T > 0, they are exponentially small in x. Apart from 
this difference, the general strategy of the asymptotic analysis is the same with however several 
interesting technical simplifications in the case T > 0. Below, we recall the general framework 
of our method. 

The 2n-fold multiple integrals of interest have the form 



\X' n 



{2my 



(27ri)" 



i=i 



A, 



• det 



1 



{A} 

{z} 



(4.1) 



r({A}) 

where Tn are holomorphic functions in some multi-dimensional strip S: |9(zj)| < a, |9(Aj)| < a, 
j = 1, . . . ,n, a being n-independent. The nearest to the real axis singularities of the integrand 



in (4.1 ) are the poles of the Fermi weight '!?(A). 



4.1 The Predholm determinant as a multiple integrals generating function 



In this subsection, we briefly recall the connection between multiple integrals of the type (4.1) 
and Fredholm determinants of integral operators [IJ. 

Let us assume for a moment that the functions J^n take the form of a pure product as 



n 



{A} 

{z} 



p=i 



(4.2) 



where </?(A) and e"^^'^^ are holomorphic in the strip |9(A)| < a. Then in (4.1), the integrals over 
the Zj separate and can be computed by taking the residues at Zj = Xj and zj = A^. A simple 
calculation leads to 



= / det[y(A,-, Afc)] d"A = d:; det[I + 7^] 



(4.3) 



7=0 
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where 



and 



27ri(A — fi) 



e 2 



(A-|.)+|{g(A)-3(M)) 



F(A) = (^(A)e-^?(^) 



(4.4) 



(4.5) 



Thus, in (|4.3|), we have identified the multiple integral (|4.1|) as the n}^ 7-derivative of the 



Fredholm determinant of the operator / + 'yV acting on M with the integral kernel (4.4). 



The large-x asymptotic behavior of the Fredholm determinant of the operator (4.4), up to 



0{e~°'^) corrections, has been established in |l44j. This asymptotic expansion as well as any of its 
finite order 7-derivatives is uniform in 7 provided 7 belongs to a sufficiently small neighborhood 
of the origin. By taking the n-th 7-derivative of this asymptotic expansion at 7 = 0, we obtain 



the asymptotic behavior of (4.3). 



In the case of interest, cannot be represented as pure product functions. However, one 
can proceed as in [l] and apply the density procedure (see |l3] for more details), that is, to 
decompose J^n on the class of pure product functions : 



n 



{A} 
{4 



k=l p=l 



(4.6) 



where the functions ipk{X) and e are holomorphic in the strip |9(A)| < a. We thus define 

Ffc(A) = ipkWe^"^^^ (4.7) 



and denote Vfc(A, fi) the kernel obtained from ^(A, fi) by the replacement F 
Then 

k=l 7=0 



Fk and g ^ gu- 



(4.8) 



We stress that in (4.8) there are no problems with permuting the symbols of integration along 



(contained in the Fredholm determinant) and summation over k induced by the density 



procedure as the convergence in (4.6) holds in the supremum norm on the strip {|9(z)| < a} 



and the weights '&{X) ensure the convergence of the integral^ 

^We would like to point out that one could avoid the manipulation of infinite sums related to the density 
procedure as follows. Namely, one first considers the case of functions J^„ that are represented as finite linear 
combinations of pure product functions and proceeds through all the steps in the re-summation below. At the 
end, once that the function one started with has been reconstructed, it is enough to observe that the answer is 
expressed in terms of a linear continuous functional defined on a dense subspace of a Banach space. Therefore, 
the obtained asymptotic expansion can be extended, by continuity, to more general classes of functions Tn. As 
this does not alter the final conclusions, we chose not to insist on that point later on. 
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4.2 Large-x asymptotic behavior of the Fredholm determinant 



We have seen that the problem of the calculation of the asymptotic behavior of the multiple 
integrals (4.1 ) boils down to the asymptotic analysis of the Fredholm determinant of the integral 
operator with the kernel (4.4). In this section, we present the results of this analysis [44]. Recall 
that the Fermi weight 'i?(A) (1.4) appearing in (4.4) is a meromorphic function in the strip 
|9(A)| < a with simple poles at A = rj", where j = 1, . . . , M. The superscripts it indicate that 
the corresponding pole lies in the upper (resp. lower) half-plane. The total number of poles 2M 
is not essential for our analysis, however we have assumed that M > 0. It is also important that 
the Fermi weight decays as a Gaussian when A — t- ±00. 

Concerning the functions g{\) and F[X) entering the kernel ( |4.4[ ), in addition to their ana- 
Ij^icity in the strip |3'(A)| < a we assume that g' {\)'d{X)F{\) and ■i?(A)F(A) are integrable along 
any curve avoiding the poles of "Q in this strip. This, in particular, means that | tr 1/| < 00. 

Now we give several definitions and introduce new notations necessary for the formulation 
of the theorem about the large-x asymptotic behavior of the Fredholm determinant det[/-|-7y]. 

In addition to the points rj", we define the points such that 1 + ^'d{s^)F{s^) = 0. For 7 
small enough, these points are slightly shifted in respect to the poles rj" (see Fig. j2|, but their 
overall number (equal to 2M) is the same in virtue of Rouche theorem. 




Figure 2: The poles of the Fermi weight rj" are depicted by o, the zeros sj" of 1 + ^'d{\)F{\) = by •. 
The integration contour ri3;i2 (dotted line) bypasses the points s\, sj, r^, rg from above and s^, S2 , 
, from below. The contour Ci3;i2 (solid line) separates tlie points Si, sj from r^, and s^, S2 
from , . 



We also introduce the sets of contours Tj-k arid Cj-k. Here j and k are multi-indices: j = 
{ji,---Jp} and k = {ki,...,kp} with 1 < js, kg < M. Also #j = #k = p, with p = 
0, 1, . . . , M. In the following, we often denote the cardinality of the sets j and k by \j\, \k\, here 
with Ij'I = \k\. 
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The contour Tj ^ is a deformation of the real axis such that, when moving from M to Tj-k, 
one only crosses the roots Sj, and the associated poles rt, r^, while all the other roots sf 
and poles rf,i^ ja or ka, are not crossed (see Fig. [2]). Note that here and below we agree upon 
the notation 



Ja -> 



{4J 



^ Ja-> 



with 



1, 



(4.9) 



Similarly, the contours Cj-k are such that when moving from M to Cj-k, one only crosses the 
roots s^, 57; all poles and other roots sf are not crossed. In particular, Cj-k separates the roots 



fit, from their associated poles rj,ri^ (see Fig. j2|. We stress that C0^0 
Finally we define an auxiliary function i^(A) as 



-0,0 



KA) 

and a contour dependent functional 



-1 
27ri 



log(l + 7^(A)F(A)), 



[ix + g'{X))u{X) dX + 



(A-/i+)2 



dX dfi 



(4.10) 



(4.11) 



where the oriented contour C is either equal to Cj-^ or Tj ^, for some j, k. In (4.11), we have 
stressed that the second integral is double. The symbol /u+ means that /i is slightly shifted to 
the left from the oriented integration contour C. 

Theorem 4.1. ^44^ Let \^\ he small enough. Then in the x — )• 00 limit and under the above 



assumptions, the Fredholm determinant of the operator I + , with V given by (4.4), admits 
the asymptotic expansion 

det[/ + 7F] = ^ exp(A:,,Jb],M)) • [l + 0(e— )] , (4.12) 
{i;fc} 

where the sum is taken with respect to all the possible choices of multi-indices j,k, including 
j = and k = %. 



Remark 1. In the following, we will use the above formula (4.12) to derive the asymptotic 
behavior of X„[J>^] from the relation (4.8 ). To do this we need to compute the n*^ 7-derivative of 



the above Fredholm determinant's asymptotic behavior (4.12). In this process one can produce 



polynomial (in the distance x) contributions to the remainders of each term in (4.12). As a 
consequence the remainders will no longer be 0(e~'^^). However, they will still be a o(e~"'^) 
for arbitrary a' < a. Hence in the following the remainders in all asymptotic formulae will be 
given up to o{e~"''^) terms for arbitrary a' < a. 



Remark 2. Moving all the contours Cj-k. in (4.12 ) to the real axis, one can recast the asymp 



totic expansion of the Fredholm determinant in the following form 



det[I + 7^] = exp(Aui[g], M) ^ Aj,k Hie 



ix{s 



b=l 



1 + e 



(4.13) 
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where Aj-j. are some constant coefficients (see |44) for tlieir explicit form). For 7 small enough, 



all the roots Sj and belong to the strip |9(A)| < a. However, certain sums ^^^li{Sj^ — s^^) 
may have their imaginary part greater then a. In such a case, the corresponding terms in the 
expansion ( 4.13[ ) (and hence, in (4.12)) can be sub-leading in respect to the correction 0(e~"^). 
Should such an effect occur, we will simply drop the corresponding contributions. 

Remark 3. In the case of interest the parameters a (and a') may depend on the coupling 
constant c only. At the same time, for 7 small enough, the roots are close to the poles of 
the Fermi weight r"^, whose positions depend not only on the coupling constant, but also on the 
temperature and the chemical potential. In particular, for h > and small T all these poles 
(and thus, the roots s''=) collapse on the real axis. Hence, for arbitrary choice of multi-indexes j 
and k one can always make ^(^[{[^^(st — )) < a provided the temperature is small enough. 



The asymptotic expansion ( |4.12 ) is uniform in 7 small enough. However, it becomes incon- 
venient for the calculation of the 7-derivatives of the Fredholm determinant at 7 = 0. Indeed, 
the contours Cj-k are squeezed between the roots s"*" and the poles r*. As — )• rj" when 7—7-0, 
strictly speaking, the contours Cj-k do not exist at 7 = 0. 

A possible way to overcome this problem is to move all the contours Cj-k in (4.12 ) to the real 
axis (see |44] ) . For our goal, it is however more convenient to deform each contour Cj-k in (4.12) 
to its associated contour Tj-k that does exist in the above limit 7—7-0. In doing so, one crosses 
the poles rt and r^. This produces additional contributions to the functional A::j-.fc([5], [j^]) and 



makes the formula (4.12) more cumbersome. 



To describe the asymptotic expansion of det[/ -|- jV] in terms of a sum over the contours 
Tj-k we introduce new notations. Let Lc[f]{X) be the Cauchy transform over a contour C of a 
function /(A): 



Lcim) 



A 



(4.14) 



Then the asymptotic expansion ( 4.12[ ) can be written in the form 



det[/ + 7^] = E ^jMial [F]) [1 + 



(4.15) 



Here, as before, the sum is taken over all the possible choices of multi-indices j ; k. The functional 
l^j-k reads 

Uj-M[9], N, [F]) = expj- lg'{XHX)dX^ fj {7'e^^'-^"'"^-^^^-^i^(^;ji^K)} • Uj,k[u] , (4.16) 



L j;fc 
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where 

Uj;kW] = { det 



Ja lib 



exp(^r.^J[0],M)) 



lil 



n ['dreg{rl)^reg{riy''^'~^^-'t'^^^'''''^i^^^^^^^^ . (4.17) 



a=l 



Above, we have introduced the residues of the Fermi weight 'd{X) at the poles rf: 'dregi'i'^) = 
Res(T?(A),A = r^). The detailed derivation of (4.16), (4.17) is presented in appendix [c| 

4.3 Asymptotic expansion of multiple integrals 

It follows from (4.8) and the from of the expansion ( 4.12[ ), that can be presented as 

^nlJ'n] = JjM^n] [l + 0(6"'^'^)] , (4.18) 

where Jj-k are some functionals that we will construct below. They can be associated with the 
contours Cj-k (or Tj-k), therefore we label them by multi-indices j,k. In order to obtain their 
explicit form one has to perform the following steps: 

• compute the n}^ 7-derivative of Uj-k{[g]j [-F]) at 7 = . 

• Substitute in this final expression F ^ Fk and g ^ Qk- 



Take the sum over k and use the decomposition (4.6) so as to express the result in terms 
of Tn or its partial derivatives. 



We stress that the derivation of this action as given in (4.19) and (4.20) is essentially a conse- 
quence of the representation (3.7) for Tn along with the reduction properties of the functions 
Wn and Vn as described in (3.16). 



Proposition 4.1. Let Tn be as given in (|3.7|) with Vn CLnd Wn satisfying to the reduction 

(4.19) 



properties (3.16). Then 



7=0 



where, for 7 small enough, 



00 ^ m „ m 

m=0 i=l i=l 



Ai>],[Vm+\j\]) ^"^L,=o- (4.20) 



Here Uj-k ([0], [v], [Vm+|j|]) is given by (4.16), where one should set g(X) = and replace the 
original functions v and F by v and Vm+|j| respectively. These new functions have the form 



v™+i,i(c.)=v„+i,i(.i 



(4.21) 
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and 



i)(a;) = Oiuj 



Tfc , {A4 

ri,{Xa + €^a}J ^-Ki 



log[l + 7^?HV„+|j|M] 



Similarly, the notation W„^_|j| means 



We remind here that rj = {rj^}^^]^^, etc. 



rt,{Aa + ea} 



(4.22) 



(4.23) 



Proof. The exphcit form of the functional Uj-if;[i'] is given by (4.17). However, to prove 
Proposition |4.1[ it is enough to know that it can be written in the following, quite general, form: 



CXD » p 



(4.24) 



where Up are some functions or distributions. Their explicit form as well as the fact that they 
depend on the poles r^^^ does not play an essential role. Substituting the above expansion into 
(4.16) and expanding the exponent of J gf'(A)z^(A) dX into a series we obtain 

m p 
d-A dP^ . Upm) Ul-g'iXaH^a)} ■ n K^a) 



°° 2|j| 
m,p=0 p 



a=l 



lil 



a=l 



X n {e'^'^^^-'^''^^^F{rl)F{r-j} . (4.25) 

a=l 



We should now differentiate (4.25 ) over 7. The n 7-derivative at 7 = of the above expression 
is non- vanishing only when n > 2|jf|, and then 



7=0 ^m,p=0{i},{g} J 



d'^A dP^ 



\j\ 



ml 



a=l 



n (-5'(Aa)^^^(A.)F^"(A.)|^j n [^'^iia)F'^^iU g^,^ 



a=l 



, (4.26) 



7=0 



where 



^0 



_1 / m p \ ^ 

— log(l + 7), Cniii}, {q}) = mV-in - 2\j\y. n n 9- ' (4.27) 



Va=l a=l 



and the symbol ^ means that the sum ov6r ^ and ^ in (|4.26l) is constrained by 
E7=iia + ELiQa = n-2\j\. 
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We have brought the result of the action of the n**^ 7-derivative into a form convenient for 
the apphcation of the density procedure. We first substitute 



(4.28) 



in (4.26). Then the sum over k is computed due to the identity 



\3\ 



k=l a=l 



a=l 
times 



a=l 
q;, times 



TT d J- { ^ ^ ^ ^ {'^ai • • • 5 a=l,...,mi \Ahi • • • ) ^b]b=l,...,p 

ii "\^+^^+^ {Aa + ea, Aa, . . . , A4a=l,...,m, {^b, ■ ■ ■ ,ib}b=l,...,p 



a=l 



. (4.29) 



ea=0 



£a — 1 times 



gf, times 



Now, we use the representation (3.7) for J^n together with the reduction properties (3.16) of the 
functions Vn and Wn- This leads to 



la times (?(, times 

J- I ^j^^k^ {Aa) • • • 1 Aa}a=l,...,m) {'^fei • • • ) '^fe}fe=l,...,p 

"^'^Z {Aa + Ca, Aa, . . • , Aa}a=l,...,m, {6> • • • , 6}fe=l,...,p 



£a — 1 times 
m 



qi, times 



b1 



6=1 



b=l 



b=l 



where the functions Vm+|j| and should be understood as in (4.21), (4.23). Substituting 

this identity into (4.26) we obtain 



m,p=0 a=l p 



m+\j\ 



lil m p 

X ni^-+iii(4)^-+biK)} n ^(^'^) n ^(^«) 



6=1 



a=l a=l 



7=0 



, (4.31) 



where 1^(0;) is given by (4.22). The sum over p in (4.31) produces lAj-k ([0], [v], \ym+\j\\) and we 
arrive at the statement of Proposition. □ 



5 Lagrange series 

We have presented the action of the functional Jj-k on the function Tn as the n*^ 7-derivative 
at 7 = of the function Gj-k{l) defined in terms of an infinite series (4.20). In this section, we 



18 



show that this series is nothing else but the continuous generahzation of a Lagrange series, in 
complete analogy with the T = case studied in |1J. It can be shown that, for 7 small enough, 
this series is absolutely convergent and that its sum can be expressed in terms of a solution of a 
non-hnear integral equation. This is a difference in respect to the T = case where the integral 
equation was, eventually, linear. 

The series (4.20) only contains the first order derivatives in respect to each €a- Once that 
these derivatives are computed, one should set all = 0. Therefore, starting from (4.20), one 
can restrict to linear in e^, a = 1, • • • , m contributions to the functions z>, Vm+\j\, and 
It means that all terms with n > 1 can be dropped., while the terms Ca eb for a^b etc should 
be kept. After such a linearization, the functions Vm_|_|j| and ^'(a') go to 

: {Aa} 



CO 



rj {Xa + ea} 



2ma - i ^ eaK{uj - Xa)^ 



and 

where 
^j-k{u} I X) 



{Aa} 



rJ {Aa + e^} 



m 

2TTia - i ^ eaK{uj - Xa)^ 



(5.1) 



(5.2) 



a=l 



X 



lil i9(ui-r+ ) 

TT ' 

a=l e ^ '■■"^ 



1 



[oj\X) 



-1 



log[l+7i9(w)$j;fc(w I X)] . (5.3) 



We recall that K{X) and 9{X) are given resp. by (1.3), (2.5). 

The linearization of the function is carried out in appendix |Bj Even though its explicit 

form is involved, this matters little at this stage of our calculations. Therefore we present it in 
a rather symbolic form 



W 



m+\j\ 



iZ + eA )-^^r,k[E'ag^JHx.y,Y.E'aebg^'\Xa,X,) 

' \a=l a=l h=\ J 



(5.4) 



where g}, (with cr = 0, it) and g^^^ are some functions. Their explicit expressions can be found 
in appendix 



B 



Here, we point out that Wj;fc depend on the three functions g^a^ {Xq)-, cr = 0, ± 
(c/ appendix [B]) . 

In all of the above notations, the subscripts indicate that the functions 0j;fc, ^j.k^ and 
Wjjfc depend on the sets of poles Vj and ry.. If these sets are empty, then we will omit the 
subscripts and write simply <I>, and W. 



Substituting (|5.1|), (|5.2|), and ^A] into (4.20) we obtain 

m 

t>j-k(>^a I 2ma - i ^ ebK{Xa 



m=0 a=l a=l 

i 7;fc 



Ab) 



b=i 



\a=l 



a=l b=l 



(5.5) 
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The above series is the continuous generahzation of the Lagrange series (c/ |45) for the scalar 
case, |46j for the multi-dimensional case, pfj and Appendix |D] for the continuous generalization). 
This series is absolutely convergent if there exists a i?o > such that for all A € Tj-j. 



1 

2^ 



log 1 + |7^(A)| 



j=i e ^ ^a ' 



(5.6) 



This clearly holds for 7 small enough and hence the series (5.5) is then absolutely convergent. 
In this case, the result of summation reads (see appendix 



2inzii, 



i)/(7^)]; 



7«(A)z,-;fc(A)dA; J g^^\X,ti)zj,kWzj,kif^)dXdf?j . (5.7) 



L j;fe 



The function Zj-k{X) appearing above is the unique solution to the non-linear integral equation 



•Zj;fc(A) — (t>j-k I 27ria — i j K{\ — uj)zj-k{^) dco^ = 0. 



(5.8) 



L j;fc 



Note that, the uniqueness and existence of this solution is provided by the convergence of the 
Lagrange series. The integral operator / -|- ICj-k acts on the contour Tj-^ with the kernel 



ICj-k{^: f^) = 'iK{X — fi)(f)'j.]^^X \ 2i7ra — i J K{X — T)zj-k{T) dr^ . 



(5.9) 



where 



(^'j.k{co\X) = dx(t>jM^\X). (5.10) 
It is easy to see that its Fredholm determinant coincides with the functional Jacobian of equation 

Observe that apart from the appearance of the Fredholm determinant described above, the 
result of the Lagrange series summation reduces to the replacement of the discrete sums involving 
the €aS by integrals over the contour Tj-k with a weight function zj-kW- Due to this fact, and 
the form of equation (5.8), one has that the original arguments (pj-k and *^j;fc of the functional 
Uj-k are replaced by zj-k and (e~^*'^^Ji'= — l)/(7'i?)) respectively (taking into account that ^j-k = 

(•g-2i7r<A,-;fe _ l)/(^t?)). 

The obtained result can be slightly simplified by deforming the integration contours Tj-k- 
Indeed, it follows straightforwardly from (5.8), (5.3) that ^-'^'^^^j.kW i^a,s simple poles at rt and 
r^. Define also the sets of zeros s+ and of e~'^'^^^^'''^^^ such that: 



g 3'"^ 3a' 



0, 



0, 



3/k 



3/k 



as 7 



0. 



(5.11) 
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We now introduce the contour Cj-k. It is a deformation of the real axis such that moving from 
M to Cj-k one only crosses the roots s'j_^, . . . , s'j^ and s^^, . . . , s^^, while the other roots sf and 
all the poles rf are not crossed (see Fig. jlj). In particular, this contour separates the poles r^/^ 
from the roots s^,, . 



Let us shift the integration contour Tj-f, to the contour Cj-k. everywhere in (5.7). Observe 
that, for 7 small enough, the pole and zero structure of the function e~'^'^'^^j;f'W in a neighborhood 



of the contours Tj-k and Cj-k. is completely analogous to the one of the function g-^^'^'^C^) KT^ 



This allows us to use the method described in appendix |C] However, now, all the calculations 
should be done in a reverse order. 



In particular, the equation (C.l) applied to the contours Fj-fc and Cj-f, and the function Zj.^, 



gives 



K{\-uj)zj.k{^)duj= J K{\-uj)zj.k{i^)duj -Y^j9{X-r^J -e{X-r+J . (5. 



a=l 



12) 



Substituting this into (5.8), we find that Zj-k{X) solves the integral equation 



-2j;fc(A) — A I 2Tria — i J K{\ — uj)zj-k{^) duj ] =0, 



(5.13) 



(recall that ((){X\X) is given by (5.3) with |j| = 0). The Jacobian of equation (5.8) naturally 

(5.14) 



turns into the Jacobian of the equation ( |5.13 ), and hence 



det[I + ICj-k] = det[/ + /C], 

^j;k Cjk 



where 



1C{X,ii) =iK{\- ij)(t)'(^\2ma-i j K{\ - T)zj.k{r) dr^ , (5.15) 



and (f)' is given by (5.10) 



Formula (C.l) also yields 



W. 



<7«(A)z,-;fc(A)dA;... ) =W( / g^^\x)zj,k{X)d\;... 



(5.16) 



(see appendix [B] for more details). Finally, making the calculations of appendix [C] in the reverse 
order we find that 



Z^i;fc([0], WA^ [(e-^---'^ - l)/(7^))]) = exp{^^. J[0], [z^.^])}. 



(5.17) 



where the functional Aa, ([0], [z]) is given by (4.11) with C = Cj-k and g{X) = 0. 
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Thus, the sum of the series (4.20) for Gj-k{'y) takes the form 



det^.J/ + /C] 



(5.18) 

We have expressed Gj-k{l) as a functional of the function Zj-kW satisfying the non-hnear 
integral equation (5.13). We now show that this last integral equation can be transformed into 
a TBA-like equation. For this, we introduce a new function Uj-k{X)' 



T 



K{\ — uj)zj-tt{uj) duo — 2'Kia + 



£(A) 
T ' 



(5.19) 



where e(A) solves the Yang- Yang equation (1.2). Then, equation (5.13) yields 

1 " 



2TTi 



log 



e(A)-u ■.^.(A) 

l + 7i?(A) e T -1 



(5.20) 



Multiplying both parts of (5.20) by iK{^ — \) and integrating with respect to A over the contour 
Cj-k, after some simple algebra, we obtain 



T 

'2^ 



Kij, - A) log 



1 + 76 T + (1 — 7)e T 



e(A) 

1 + e T 



dX. (5.21) 



It is easy to see that the roots and Sf, satisfy the equations 



l + 7exp(-n,-;fc(s±/J/T) +(l-7)exp(-e(s±/J/r) =0 



(5.22) 



When deforming the contour Cj-k to the real axis in ( 5.21^ we cross these roots, but do not cross 
the poles of the Fermi weight rf, i.e. the points where 1 + exp(— e(r^)/T) = 0. Therefore, we 
obtain 



Uj-k{lj) - e(/^) + ^TiiaT = j K{p, - A) log 



"3;fc(^) ^(A) 

1+76 T -|-(i_^)e T 



d\ 



+ — / K{fi-X)log 



1 + 6 T 



dX-iTY, [0{sl -lA- 0{sl^ - IJ)] . (5.23) 

a=l 



Finally, using the Yang- Yang equation we arrive at 

T 



Uj-kifJ-) = — (/i + l-KiaT) 



2tt 



i^(/i-A)log 



1+76 T + (1 — 7)6 T 



dX 



iT 



n 
a=l 



^i)-9m^-fi) . (5.24) 
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Thus, we have reduced equation (5.13) for the function Zj.fc(A) to equation (5.24) for the function 



Uj-k{X). In the following, we shall consider the function Uj-kW as the primary object, while the 



function Zj-fc(A) will be defined by (5.20). 

Let us summarize the obtained results. We have expressed the large x asymptotic behavior 
of the multiple integrals In[J~n] (4.1) in terms of the function Gj-kil) via (4.18), (4.19). The 



function Gj-kij) is given by (5.18). It is a functional of Uj-kW, the solution to the non-linear 



integral equation (5.24). In its turn, Uj-k{^) depends on the sets of parameters and 



satisfying to the condition (5.22) 



Observe that the roots st/i. depend on 7: 



deformations of the poles rjy^^ such that 5^/^(0) 



s^lkil)- treat them as 7- 



3/k- 



Similarly, the function Uj-k{^) 



Uj-ki^i 7) can be considered as a 7-deformation of the function Uj-k{\^ 0) which satisfies equation 
( [fl4| at 7 = 0: 



nj.fc(A,0) = e(A) - 2maT 



a=l 



A) - e{r,^, - A) 



(5.25) 



The above 7-deformations are analytic in 7, at least for 7 small enough. 

To conclude this section we simplify the notations used above. Recall that we have originally 
used the multi-indices j and k in order to denote certain subsets of the Fermi weight's poles rj" 
and r^. Let us enumerate all subsets of multi-indices {j; h} by one number, say i. Then every 



subset T^ik can be enumerated as In particular, we agree upon {r"'=}c 



However, in all 



other respects, the order of enumeration is not essential. Given a subset {r*}j we can uniquely 
determine the roots and the corresponding function Ui{\) as analytical 7-deformations of 



{r*}j and Ui{\, 0) (5.25). In its turn, given we can define the integration contours Cj, and 



hence, find the function Gi{j). Thus, the asymptotic behavior of the multiple integrals (4.1) 
can be written in the form 



1 + e~ 



7=0 



(5.26) 



where i enumerates the subsets of multi-indices j and k. 



We do not write explicitly the upper limit of summation over i in the formula (5.26). It 



depends on how many poles of the Fermi weight belong to the strip |9(A)| < a'. In particular. 



it follows from the Remark 3 given on page 15 that for h > and T small enough and arbitrary 
io there exists T(zo) such that 



lim e'^"a"Gi„(7) 



7=0 



0, for T < T{io). 



(5.27) 



Therefore in this case the sum in (5.26) may contain an arbitrary number of terms. 
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6 Asymptotic behavior of the correlation function 



We have presented the large x asymptotic behavior of the integrals Xn[J-n] (4.1) in the form 



(5.26). We should now substitute this result into the series (3.6), what leads us to 



oo ^ oo ^ 

n=0 - — " ■ 



n=0 



7=0 



(6.1) 



Since Gi{'y) does not depend on n, the first series gives the Taylor expansion for Gi{'y) at 7 = 1. 
Hence, this series should result in Gj(l) provided it is convergent at 7 = 1. One can easily 
convince himself that it is true at least for the case c = 00. Indeed, in this case the sum of Gi(l) 



coincides with the asymptotic expansion of the Fredholm determinant (3.15) obtained in |44) . 
Therefore considering the QNLS model with c < 00 as a smooth deformation of the free fermion 
case we assume that this series is convergent at 7 = 1, and thus we obtain 



1 



Z^Z^ ,-9^^,(7) 



n=0 



7=0 



E 



G^{l) 



(6.2) 



Note that the convergence of the series (6.2) is related to the convergence of the Lagrange series 



(5.5) at 7 = 1. This means that the analytic 7-deformations described in the end of the previous 
section can be continued from a vicinity of 7 = to the point 7 = 1. Hence, for the computation 
Gi{l) it is enough to set 7 = 1 in the equations for the functions Ui{X) and Zi{X). The non-linear 



integral equation (5.24) for nj(A) takes the form (2.4), while the equation ( |5.20 ) for Zi{X) turns 
into 

1 , / 1 + e" 



T 



l + e 



£(A) 

T 



Substituting this into the functional ^(^.([0], [zi]) we find 
^c^iM, [zi]) = -ix / z^{X) dX + 



(A-M+)^ 



(6.3) 



(6.4) 



Ci 



and hence, every Gi{l) is proportional to e~^^' with pi given by (2.9). Note that the generating 
function (e^'^*"^'')T' is not necessary a decreasing function of x at arbitrary complex value of a. 
Therefore the real parts of the obtained pi are not necessary non-negative. However we should 
have ^{pi) > at Q = 0, since this case describes the physical two-point correlation function 
{j{x)j{0))T- This property is confirmed by numerical computations (c/ appendix [A|) , and the 
analysis of the low-temperature limit |47) . 



Concerning the series corresponding to the reminder in (6.1), we assume that it remains 



exponentially small with respect to the first term, provided the sum over i in (6.2) is restricted 



by some zq. Similarly to (5.26) we do not write this upper limit of summation explicitly, because 



we can not find its exact value for arbitrary temperature. One can show however (c/ jlT]) that 
io grows as the temperature decreases. In particular, it goes to infinity when T — )• 0. 
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Thus, we reproduce the expansion ( |2.8[ ) 



X — )• oo. 



(6.5) 



where Pr, 



maxi(JR.{pi)). 



It remains to describe the constant coefficients B[ui] in (2.8). One part of these coefficients 
comes from the constant term in the functional ([0], [zi]) (6.4). Another part is equal to the 
Fredholm determinant det^ [I + /C] appearing in equation (5.18). Observe that, at 7 = 1, the 



derivative (p' (5.10) reads 



e(A) - Ui{X) 



2717 - "i(^) 

1 + e — ~ 



Thus, 



det[/ + /C] = det 



2tt 



where the kernel is the analog of K^'^^ where e has been replaced by 



Uj-. 



"i(M) 

1 + e T 



(6.6) 
(6.7) 

(6.8) 



This operator acts on the contour Cj and not on M as it was the case for / + K^^'' . 

Finally, the most complicated part of the coefficients B[ui\ comes from the function W (this 
part of the computations is given in appendix p|). Combining all these results we get 



1)^ exp 



B[ui 



Zi{X)zi{fi) 
;(A-^+)2 



dX dfi — Co[zi, Ci 



r L^-^ [zi]{8i+ic) ^Tvia+L^^ [zi]{d\~ic)^ ,~L^^ [zi]{92-ic) ^2TTia-L^_ [zi]{02+ic)-. 



det 



I+^^U^%^] det /+ J^C/(2)[^,] 



det [I - ^K(^)] det [l - ^Ki^^ 



(6.9) 



Recall that L^ [zi] is the Cauchy transform (4.14) of the function Zi{X) on the contour Cj. The 



functional Cq as well as the kernels of the integral operators {w , w' , [zi]) and U^'^\w,w',[zi]) 



are given in appendix [B] (see resp.(B.6), (B.9), (B.IO)). The latter operators act on an counter- 



clockwise oriented contour surrounding Cj. 



7 The leading term and corrections 

Taking the second a-derivative of the expansion 

for (e2™2.^^ at a = we obtain the 
asymptotic behavior of the density-density temperature correlation function (j(x)j(0))T- We 
show below that the leading term of this asymptotic expansion is produced by e~^'^° correspond- 
ing to the contour Cq = M while the other pj's lead to sub-leading corrections. 



25 



7.1 The leading term 

As expected on general grounds, we obtain the following leading term of the asymptotic behavior. 
Proposition 7.1. The leading term of the asymptotic behavior of the two-point function 

(j-(^)i(0))T = {jix)jm)? + 0{x-^) . (7.1) 
is given by the square of the density of particles in the state of thermal equilibrium: 

(j(x)j(O))?) = ( 1 p,{X)dXy = (7.2) 
This contribution stems from the choice of the contour Co = M. 



Proof. Setting Co = M in (2.7) we obtain 



no(^) = ^^-ho^-— K{ix-\) log (l + e-— j dX. (7.3) 

R 

The form of this equation coincides with the one of the Yang- Yang equation up to the shift of 
the chemical potential. Thus, we have 

uo{\) = e{\\had, (7.4) 



where we have insisted explicitly on the dependence of e{\) on h. Substituting this into (2.9) 
we find 

1 f /l + e"=(^l''")/^\ 1 



^° = y l + e-<m/T )dX=^{V-V^), (7.5) 



where we have defined 



7. = ^ y log (l + e--(^l'^)/^) d\ . (7.6) 



The constant Va in (7.5) corresponds to the shifted chemical potential, i.e. it is defined by (7.6) 
with e{\\h) replaced by e{\\ha). This result reproduces the prediction of |12| . The quantity V 
can be interpreted as the pressure in the gas, whereas Va would correspond to the pressure in 
the presence of a shifted complex valued chemical potential. 

Thus, we see that po ~^ as a — ^ meaning that the contribution stemming from the 
choice z = (that is [j; k) = (0,0)) to the density-density correlation function does not have an 
exponential decay. Applying the second order derivatives over a and x to the presumed leading 
order term of ^e^'^*"®'')^ we obtain 



(7.7) 

a=0 
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and thus, to proceed further we should find uq{X) and 2;o(A) up to hnear in a terms. Differenti- 

(7.8) 



ating (7.3) over a at a = we obtain 



da 



a=0 



duo{X) 



dX, 



a=0 



and hence, 



'Uo(A) = e(A) - 4TT^iaTpt{X) + ©(a^), 



(7.9) 



where pt{X) is the total density satisfying the integral equation (1.5). Equation (1.5) can be 



solved in terms of the resolvent R{X, fi) defined by (/ — ^K^^^^ (I + R) = I: 



2Trpt{p) = l + y ^(/^>A) dX. 



(7.10) 



Substituting (7.9) into (6.3) we obtain 

zoiX) = -27rapt(A)i9(A) + ©(a^) 



and thus, due to (2.9), 



dpo 
da 



-27ri 



a=0 



Hence, 



J pt{X)^{X) dX = -2^i J pp{X) dX. 



Pp{X)dX • B[uo]|„=o, 



(7.11) 



(7.12) 



and it remains to apply Lemma 7.1 in order to conclude that B[uq] = 1 at a = 0. 

Lemma 7.1. B[uq] = 1 at a = 0. 

Proof. Since zq^X) — )• at a — )• 0, we conclude that 

expf / ^o(^)^o(^) dXdf, - Co[zo,R]] = 1, at a = 0. 
\Jr (A - p+r J 



(7.13) 



□ 



(7.14) 



For the calculation of the limit of the remaining part of equation (6.9), we use the obvious 
properties of the Cauchy transform L]g[zo](a;): 

Lk[zo]{u + iO) - Lk[zo]{uj - iO) = 2TTizo{uj), 

L]^[zo]{u} + ic) - Lik[zo](w - ic) = ^(no(A) - e{X) + 2-niaT) = 27ria(l - 2Tipt{X)) + O(a^). 

(7.15) 
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The last property follows from (5.19). Using (7.15) we find at a — )• 

g27ria 



and 



g(i*o(ei)-£(ei))/T _ ^ 
I „2ma ^-2ma+L^[zo\(e2-ic) 



(7.16) 



(7.17) 



Finally we should compute the ratio of determinants in the second line of (6.9). Consider, 
for example the determinant of the operator /+ ^U^^\ Due to the factor e^*[^'']("'\ the kernel 
U^^\w,w') has a cut on the real axis. Hence the action of the integral operator can be 
reduced to an action on the real axis with the kernel 



[>(i) 



r(R) 



f/d) 



(7.18) 



where u!^^ are the limiting values of ?7^^^ from the upper (lower) half-planes. Using equat: 
(7.15 ) we obtain 



;ions 



det 
r(K) 



det 



(7,19) 



where the operator in the r.h.s. acts on M and its kernel reads 



{uo{w)-eiw)) /T _ ^ 



Setting now a = and using (7.9), (7.11) we obtain 

= det 



det 



1 



«=o 



(7.20) 
(7.21) 



Thus, in terms of the resolvent i?(A, p) the ratio of determinants can be presented as 
detail - i^{K('Hw,w') - K(^-\ei,w'))] 



detail -^K(^){w,w')] 



det[I + R{ei,w')] 



1 + j R{ei,w)dw = 2TTpti9i), (7.22) 



where we have used (7.10). Similarly one has 

detail -^K(-){w,w')] 
and we arrive at the statement of Lemma 17.11 



2vrpt(^2), 



a=0 



(7.23) 
□ 
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7.2 The corrections 



It is clear that for other terms in the expansion (2.8) Ui{X) ^ e{X) even at a = 0. This yields 
Pi ^ at a = and C 7^ M. Therefore after taking the second a-derivative and setting a = 
the corresponding terms of the expansion contain exponential factors e~^'P\ Our numerical 
computations together with the low-temperature limit support that the real part of the pi's 
are indeed positive and, hence, the corresponding terms in the asymptotic expansion are sub- 



leading. In particular, in the low-temperature limit, equation (2.7) can be solved for all possible 
choices of subsets j,k. Respectively, one can find all pi. It follows from this analysis, that 
the main sub-leading contribution to the asymptotic behavior of <|e^'^*°®^)T originates from the 
contours Cj, where {s"*"}! are 7-deformations of the poles {r*}, lying nearest to the real axis. 
We conjecture that this property holds for the finite temperature as well. 

Finally, we would like to give an interpretation of the results obtained from the viewpoint 
of the quantum transfer matrix approach [SS]- As we have mentioned already this method was 
developed for the description of quantum spin chains, but it was shown in |32| that in a special 
continuous limit it can be applied to the description of thermodynamics in the one-dimensional 
Bose gas. The central object of this method is a function ao(A) constructed via certain solution 
of Bethe equations and satisfying a non-linear integral equation. Knowing this function allows 
one to calculate the maximal eigenvalue Aq of the quantum transfer matrix Tq. In the continuous 
limit mentioned above, one can establish the following correspondence 

ao(A) ^ e-W/^, 

1 r. r -sW/T\^^_ (7.24) 



Ao^^yiog^l + e ^ 

R 

The sub-leading eigenvalues Aj of Tq can be expressed in terms of functions aj(A) satisfying 
the same type of the integral equation, but on deformed integration contours surrounding some 
zeros of 1 -|- ai(A). It is easy to show that, in the continuous limit, these contours go into the 
contours Cj. Then one should have the following correspondence 

^ i_ I log (1 + e-«»(^)/^) dX. (7-25) 



Taking into account (2.9) we obtain the following mapping 



Ao{h) 



^e-^f% (7.26) 



where m is the lattice distance that scales to the distance x in the continuous limit. In (7.26), we 
have also stressed that the eigenvalue Ao(/i) corresponds to the quantum transfer matrix Tq{h) 
with the chemical potential /i, while the eigenvalues Ki{ha) correspond to the quantum transfer 
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matrix Tq{ha) with the shifted chemical potential ha- Thus, the rates of the exponential decay 
in the asymptotic behavior of (e"®'')^ appear to be nothing else but the ratios of the eigenvalues 
of the quantum transfer matrix (in the continuous limit). 



Conclusion 

The main goal of this article was to derive the long distance asymptotic behavior of temperature 
correlation functions of one-dimensional bosons in the framework of the algebraic Bethe ansatz 
approach. Considering the example of the two-point function of densities we showed that this 
asymptotic behavior can be expressed in terms of solutions of the TBA equations. This result 
is quite expected from the viewpoint of the QTM method. 

Comparing to the QTM approach, we should say that our derivation is certainly more in- 
volved technically. However in our opinion these two methods are complementary, in particular 
speaking about their range of applications. Our starting point is the master equation, which 
is a multiple integral representation for various correlation functions, including the dynamical 
ones. Recall that such type of representations is known now for an arbitrary integrable model 
possessing the six- vertex i?-matrix |6]. The expansion of the master equation into the series 
allows us to take the thermodynamic limit. Then in the asymptotic regime this series turns 
into the generalized Lagrange series, which can be summed up explicitly in terms of solutions to 
some integral equations. The last ones appear to be just the equations providing the dressing 
of bare variables originally entering the representation for the correlation function. We have 
demonstrated such a mechanism already in the case of zero temperature correlations [1]. We 
have now shown how it works in the case of thermal correlation functions. 

Our derivation however is not free of several assumptions of technical kind. For example, such 



questions as the convergence of the series (3.6) and (6.2) formally remain unsolved (although 



these series are absolutely convergent at the free fermion point). We think nevertheless that 
the results obtained under these assumptions are themselves strong arguments in favor of their 
validity. As an additional test for our results we will consider the low-temperature limit of 
our asymptotic expansion in a forthcoming publication. In this limit, the non-linear integral 
equations determining the correlation functions in the asymptotic regime can be solved explicitly 
in terms of dressed physical quantities such as the energy, momentum, and charge. In particular, 
we will show |47] that the low-temperature limit of our results reproduces the conformal field 
theory predictions for the exponential decay of the correlation functions. We are however able 
to go well beyond these predictions as our results not only hold for all ranges of temperatures 
but also provide explicit formulae for the corresponding amplitudes. Moreover, it will be shown 
that the dependence of these amplitudes in terms of powers of the temperature in the low- 
temperature limit can be computed explicitly, in complete agreement with the CFT predictions 

m 
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A Numerical analysis of the correlation lengths 



The inverse correlation lengths ^{pi) are defined in terms of the integral (2.9) involving the 



solution Ui to the TEA non-linear integral equation (2.7). Hence, the pi depend on the coupling 



constant c, the temperature T, and the chemical potential h. Apart from these parameters, 
they also depend on the choice of the system of roots and it is to this dependence that 

the subscript i refers to. In this appendix, we gather four plots resulting from our numeri- 
cal computations of several correlation lengths (the real part of the pi) as functions of their 
parameters. 




h 




Figure 3: The dependence p = p{h) in a) for a fixed temperature T — 2 and different coupling constants 
c — 7 (black), c — 10 (blue) and c ~ oo (red). In b) the coupling constant is fixed at c = 10 and the 
temperature is set to T = 1 (black) . T = 3 (blue) and T — 5 (red) . 

The situation considered in the first three figures (Fig. |3]a) and b). Fig. |4]a)) corresponds 
to the correlation length determined in terms of the solution to the non-linear integral equation 
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a) b) 




Figure 4: In a), the dependence p = P{T) is presented for the coupling constant c = 10 and at 
different values of the chemical potential h = —\ (black), h ~ A (blue) and h = 9 (red). In b) we have 
represented ^{p) as a function of the chemical potential h for the coupling constant fixed to c = 10 and 
the temperature to T = 2. The three cases correspond to one pair of roots in various positions : (i) 
is the nearest root in the first quadrant, s~ is the nearest root in the fourth quadrant (black); (ii) s+ 
is the nearest root in the first quadrant, s~ is the nearest root in the third quadrant (blue); (iii) is the 
next after the nearest root in the first quadrant, is the nearest root in the fourth quadrant (red). Two 
cases with two pairs of roots s*: (iv) sf are the nearest roots in the first and the fourth quadrants, 
are the next after the nearest roots in the first and the fourth quadrants (violet); (v) sf are the nearest 
roots in the first and the fourth quadrants, §2 are the next after the nearest roots in the first and the 
third quadrants (brown). 



with one root in the first quadrant of C and one root s~ in the fourth quadrant. These roots 
are built as 7 deformations of the poles of the Fermi weight lying closest to the real axis and 
located in these quadrants. For this case, the corresponding p is real. 

On Fig. [I a) we have represented p as a function of chemical potential for a fixed temperature 
and several values of the coupling constant c. Similarly, in Fig. [3] b) we have also represented p 
as a function of chemical potential but for a fixed value of the coupling constant c and different 
values of the temperature. On Fig. [l] a) we have represented p as a function of temperature for 
a fixed value of the coupling constant and several values of the chemical potential. For all these 
quantities, the numerical values are given in dimensionless units after proper rescaling of the 
TBA non-linear integral equation (2.7). 

The last plot Fig. |4] b) represents different correlation lengths ^{p) as a function of the 
chemical potential and for a fixed value of the coupling constant and the temperature. There, 
we have considered several possible configurations of the system of roots {s^} for which p are 
complex numbers. 
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B Linearized form of Wj, 



The explicit expression for the function Wm is given by (3.10)-(3.13). It contains Fredholm 
determinants whose kernels depend on the products (3.8). At {zk} = {Afc + e^} we have, to the 
linear order in each of the e„'s 



lil 



w — r, + ica 



w — rj + tea 

a=l Ja 



exp £ 6.5«(AJ 



a = 0,±, 



(B.l) 



.a=l 



where 



1 

w — X + iac 



0- = 0,±. 



(B.2) 



Apart form the products of the type (B.l ), the function Wm also contains the product of the 
functions y^"^^'^'^ (see (3.10)). Setting {zjS\ = {A^ + e^} we obtain 



n 



n 



{m |j| 
a=l b=l 

X n !:^::^:Sl:j::5:S -e^p(EE---^^"^(A.>A4, (b.3) 



=1 -^fc, 



,a=l 6=1 



where 



{X — fj, — ic) 



2 ■ 



(B.4) 



The expressions (B.l), (B.3) should be substituted into (3.10), what gives us the linearized 
form Wj-k ( |5^ . 

The summation of the Lagrange series leads to the replacement of the sums by integrals over 
the original contour Tj-k 

m „ 

] — 1 ^ 



k=l 

m m 



E E £fc^fc'5'^^^(Afc,Afc') ^ / g^'^\X, fi)zj.kiX)zj-kifJ-) d^dfi = -Co[zj.k,'rj-k], 



(B.5) 



k=i k'=i 



where Zj.k{X) solves the integral equation (5.8), -f^r^-.j, is the operator of Cauchy transform (4.14) 
and for any contour C 

Zj-kWzj-kil^) 



Co[zj-k, -C] 



{X — fi — icy 



dXdfi . 



B.6) 
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Substituting these formulae into (3.10) and moving the contours from Tj k to Cj-k we get 



that all the explicit dependence on the poles r^y^ cancels out and 

(A) 



W / g^J\X)zj,kWdX; g^^\X, f,)zj.kiX)zM dXdfi) , (B.7) 



where 



W N / g^^\x)zj,u{X)du \ ■ g^''\X,^l)zj,k{\)zr,k{^J)dXd|Ji 

e-Co[.,;fe,C,;,](g2.ia _ (/ + 2^f/(i)[2,-.fc]) det (/ + ^fi^^\zj,k\ 



fB. 



r La [zj-k\{di+ic) 2TTia+LA [zj.k](0i-ic)^. -La [zj.k\{62-ic) 2-Kia-LA [zj-k\{92+ic)^ 
g ''3;fe •" — g ''j;fe g ''j;fc — g ''3;fe 

and the kernels of the integral operators U^^\w^w' ^ [zj;k\) and JJ^'^\w,w' , [zj-k\) have the form 

g ''j;fe — g '-Jife 

g ''j;fc — g 

Both integral operators act on a counterclockwise oriented closed contour surrounding the con- 
tour Cj-k- 

C Deformation of the contours Cj k 



In this section we calculate the difference Ac^.i^Hg], [i^]) — -^rj.kiidjy W]) (see (4.11)). We shall 
consider some fixed contours Cj-k and Tj-k- Hence, for brevity, we omit the subscripts j;k in 
the following. Moreover, it is clear that, without loss of any generality, we can set j = {1, . . . ,n) 
and k = (1, . . . ,n). Also we note that we can carry out the intermediate computations up to 
integer multiples of 2iTT. This is justified in as much as we take the exponential at the end. 

Observe that the contour C U —T surrounds the points in the counterclockwise direction 
and the points in the clockwise direction. Therefore if /(w) is holomorphic in a domain 
containing T and C, then 

I nu)v{u)du = j f{uMu;)du; + j2{firk)-fi4)) ■ (d) 
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Using (C.l) one can easily calculate the difference of single integrals entering (4.11): 

- j\ix + g'{\))v{\)d\ + j\ix + g'{\))y{\)d\ = Y,{i<rl-Tl)+g{rl)-g{rl)). (C.2) 



The calculation of the difference of the double integrals in (4.11) is more involved. It is 
convenient to present i^(A) in the form 



1 



fr(A-4)(A-.,-) 



^,Vi(A-r,+)(A-r,) 



where i'(A) is holomorphic in a domain containing both contours C and V. Let 



^ A ^ 

— -^dXdfx, J{T) 

{X- fi+y 



(A-/z+)2 



dX dfj,. 



(C.3) 



(C.4) 



Let us also introduce auxiliary functions 



a- 



-(A;C) = log|n^ 



k=l 



a. 



-(A;C) = log|n^ 



k=l 



+ I ' 



and 



(C.5) 



(C.6) 



It is easy to see that the functions a±(A; C) are holomorphic to the left (resp. to the right) from 
the contour C. Similarly the functions a-|-(A;r) are holomorphic to the left (resp. to the right) 
from the contour F. All the functions (C.5) and (C.6) behave as 0(A^^) when A — )• oo, and 



f " (X — s^)(X — s^) 

a+(A; C) + «_ (A; C) = a+(A; T) + «_ (A; T) = log m )t ^ 

1 (A-r+)(A-r. ) 



(C.7) 



Using (C.3) and (C.5) we can present J(C) as a sum of four integrals J{C) = Jj{C), 



where 



Ji(C) 
Ja{C) 



P(A)P(^) 



dX dfi, 



1 

'2TTi 

1 

'2TTi 



(A-^+)2 

t>(A)(a+(/x;C) + a_(^;C)) 
(A-^+)2 

(a+(A;C) + a-{X;C))i'{fi) 



dX dn, 
dX djj,, 



1 

(27ri)2 



(A-M+)2 

(a+(A;C) + a_(A;C)) (a+(/u;C) + a_(^;C)) 
(A-A^+)2 



dX dfi. 



(C.8) 
(C.9) 
(C.IO) 
(C.ll) 
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Similarly J(T) = Yl^=i "^jO^) where Jj{T) is obtained from Jj{C) by replacing everywhere C by 

r in ([a8l)-( |cTTT| ). 

Since i'(A) is holomorphic in a domain containing the contours C and F, we conclude that 
Ji(C) — Ji(r) = 0. The integrals J^iC) and J4(T) can be taken explicitly. For example, using 
the analytic properties of a±{X;C) we have 



1 



1 



a 



+ (At; C) («+(//; C) +a_(/x;C)) d/x = — / a'^{ij;C)a-{ij;C) dfi 



-f±(- 



a_(/u;C) d/i = log ]^ 



. j,k=l 



. (C.12) 



The expression for Ja{T) is obtained from (C.12) via the replacement r+ -H- r . Then we have 

(C.13) 



j4(c)- j4(F)=iog<j n 



Calculating J2 and J3 we can take explicitly only one of two integrals. We have 

1 



J2(C) = I a'_(A;C)P(A)dA = / E (a^ 



^k ^ - ^fc 



^ z>(A)dA, 



(C.14) 



and 



J3{C) 



^k f^-H 



D{^i)dn. (C.15) 



The expressions for J2(F) and J3(F) can be obtained from (C.14) and (C.15) via the replacements 
C — > F and r"*" -H- r~. Then we have 



1 



^k ^ - Sfc 



P(A) dA 



cur '==1 



^ z>(A)dA. (C.16) 



The integral in the first line of (C.16) vanishes as the contour C U — F does not surround the 
points s^. The second integral in (C.16) gives 



s=2,3 p fc=l ^ 



^k 



>(A) dX - ^(z>(r+) + D{r,)). (C.17) 



k=l 
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Here we have used that the contour C U —T surrounds the points rt in the counterclockwise 



direction and the points in the clockwise direction. Substituting into (C.17) the function 



t'(A) in terms of z^(A) via (C.3) we obtain after simple algebra 

n 

{JsiQ - Js{T)) = Y,{2LrW]{r,) - 2LrMK) " 2vrii>(r+) - 27ri9{r~,)) 



s=2,3 



k=l 



+ 21ogn l-^-'^], (C.18) 



where we have used the Cauchy transform over the contour T of the function iy{X) (see (4.14)). 



It remains to find t'(?^^). These numbers can be expressed in terms of the residues of the 
Fermi weight i^regi^^) points r^. We have 



Hence, 



1 + 7^(A)F(A) = e---W = e---W ^l^l^^ 



1 



(C.19) 



(C.20) 



Thus, combining (C.13), (C.18), and (C.20) we obtain 

2 



,m-JiT)^ det 



1 \ " 

1 k / k=l 



where we have used 



det 



j,k=i 



(C.21) 



(C.22) 



Taking into account (C.2) we immediately arrive at the formulae (4.16), (4.17). 



We insist that all the above computations have been done by using the sole pole/zero struc- 
ture of e-2i7n/(w) 'pj^gy ^j.g thus valid in the specific case where one considers v = zj-k. 

D Continuous generalization of the multiple Lagrange series 

The continuous generalization of the multiple Lagrange series has the form 



n=0 j j=l j=l \a=l J \a=l J 



fD.ll 



e,=0 
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Here the integrals are taken over some contour £, the functions /, F, ^, and hi are holomorphic 
the corresponding neighborhoods. The argument of the function F may also contain double 



sums 5 /i2(Aa, Aft)eaefe etc, like, for example, in (5.4). Apart from obvious modifications, 
the presence of such multiple sums does not affect the result. Hence, we have omitted these 
arguments of the function F for brevity. 

This series was studied in [1] . If there exists a Rq > such that 



sup sup 

(/)e[0;27r] IJ.&C 



fD.2) 



then it is absolutely convergent. In such a case, the result of the summation reads 
G= 



det£ 



where the function z{^) is the unique solution to the integral equation 



2(/^) = /( j e(A,M)^(A)dA 



(D.3) 



(D.4) 



The denominator of (D.3) contains the Fredholm determinant of a linear integral operator acting 



on C. If the function F depends on multiple sums, then these should be replaced by multiple 
integrals, like, for example. 



n „ 



(D.5) 



a,b=l 



etc. It is important to note that the form of the integral equation (D.4) does not depend on the 
function F. 



The result (D.3) can be directly applied to the summation of the series (5.5), where the 



integration contour £, coincides with Tj-k and the function / is (pj-k (5.3) 
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